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We numerically study disorder effects in Bernevig-Hughes-Zhang (BHZ) model, and find that 
Anderson transition of quantum spin-Hall insulator (QSHI) is determined by model parameters. 

The BHZ Hamiltonian is equivalent to two decoupled spin blocks that belong to the unitary class. 

In contrast to the common belief that a two-dimensional unitary system scales to an insulator except 
at certain critical points, we find, through calculations scaling properties of the localization length, 
level statistics, and participation ratio, that a possible exotic metallic phase emerges between a QSHI 
and a normal insulator phases in InAs/GaSb-type BHZ model. On the other hand, direct transition 
from a QSHI to a normal insulator is found in HgTe/CdTe-type BHZ model. Furthermore, we show 
that the metallic phase originates from the Berry phase and can survive both inside and outside the 
gap. 

PACS numbers: 72.15.Rn, 73.20.Fz, 73.21.-b, 73.43.-f 


I. INTRODUCTION 

Topological insulators (TI), which are identified as a 
class of quantum state of matter, have generated in¬ 
tensive interests recently^— The two dimensional (2D) 
TI, i.e., quantum spin-Hall insulator (QSHI), is charac¬ 
terized by odd pairs of counter propagate gapless edge 
states.—^ The QSHI was firstly realized in HgTe/CdTe 
quantum well (QW)^£ and subsequently in InAs/GaSb 
QW£Z These two experimentally realized QSHI systems 
can both be represented by the inverted bands Bernevig- 
Hughes-Zhang (BHZ) model but with different param¬ 
eters. Specifically, the coupling strength between two 
inverted bands in InAs/GaSb QW is an order of magni¬ 
tude smaller than that of HgTe/CdTe QW. Considering 
the remarkable parameter difference, one open question 
remains, namely whether this difference will have some 
physical consequences in these two QSHI systems. 

The Anderson metal-insulator transition in 2D dis¬ 
ordered system manifests as a lasting research issue in 
condensed-matter physicsA-ii Generally, the disordered 
electron systems can be classified into three universality 
ensembles according to the random matrix theor y 14 i 15 
In the presence of time-reversal symmetry (TRS), the 
system is classified as an orthogonal ensemble if spin ro¬ 
tation symmetry is preserved; otherwise, it belongs to a 
symplectic ensemble. In contrast, when the TRS is bro¬ 
ken, the system turns into a unitary ensemble : 14 ' 15 For 
a TRS QSHI with Rashba spin orbital coupling (SOC), 
the system falls into the symplectic ensemble because 
spin rotation symmetry is broken. In contrast, without 
Rashba SOC, the QSHI is divided into two spin species 
of quantum anomalous Hall (QAH) systems, which be¬ 
long to the unitary ensemble. Previous studies on metal- 
insulator transition in a QSHI can be summerized into 
two paradigms: (i) for (symplectic) QSHI with Rashba 


SOC, it was found that TI and normal insulator (NI) 
phases are separated by a metallic phased - — (ii) for 
a (unitary) QSHI without Rashba SOC, a direct tran¬ 
sition from a TI to a NI was discoveredSurpris¬ 
ingly, recently a crossover from weak localization to weak 
anti-localization (WAL) is suggested in the BHZ model 
without Rashba SOC<22r— Since WAL can add a positive 
correction to the /3 function^ ' 23 ' 24 there is the possibility 
that a metallic phase might exist in the 2D unitary sys¬ 
tem (QSHI) with weak disorder, although that is against 
the traditional view. 

In this paper, we study whether a metallic phase can 
exist between TI and NI phases in a unitary QSHI sys¬ 
tem. Starting from the BHZ model Hamiltonian^* we 
calculate the localization length and two-terminal con¬ 
ductance numerically, and we analyze the scaling behav¬ 
ior of the system : 10 ' 25 It is worth noting that the BHZ 
Hamiltonian is equivalent to two decoupled spin blocks 
which belong to the unitary class. The main results are 
summarized in the phase diagrams in Fig. [l] In all these 
phase diagrams, we find a metallic phase between TI and 
NI phases, in contrast to the common view of Anderson 
transition behavior in a 2D unitary class. Furthermore, 
we find that different parameters in BHZ model (uni¬ 
tary system) can lead to different Anderson transition 
behaviors. The transition from a TI to a metal is likely 
to exist in InAs/GaSb-type BHZ model but not in the 
HgTe/CdTe-type BHZ model [see dash lines in Fig.[ljb)]. 
By employing the Berry phase, the parameter-dependent 
metallic phase can be well explained: 22 ' 26 

The rest of the paper is organized as follows. In Sec. 
II, we introduce the lattice model Hamiltonian and give 
the details of numerical simulations. In Sec. Ill, we show 
the main results by scaling of the localization length and 
the conductance. We discuss energy level statistics and 
participation ratios of eigenstates in Sec. IV. In Sec.V, 
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(a) A=0.28, Ep=0 


(b) M=0.38, Ep=0 


(c) A=0.28, M=0.38 



FIG. 1: (Color online). The phase diagrams for disorder strength W and (a) topological mass M, (b) electron-hole hybrid 
strength A, and (c) Fermi energy Ef■ The black filled circles are critical points determined by finite size scaling, while the 
red filled squares for W = 0 are analytic critical points. In (b), the dark cyan dash lines (InAs/GaSb-type BHZ parameter) 
indicates the TI-metal-NI transition. The dark yellow dash lines (HgTe/CdTe-type parameter) shows direct NI-TI transition. 


we discuss the phase diagram and interpret our numerical 
results by the Berry phase. In Sec.VI, a brief summary 
is presented. 


II. MODEL AND METHODS 


We consider the disorder BHZ Hamiltonian on a square 
lattices^^ 


H = J2^Ei^+ J2 v\T a <p i+a + H.C., ( 1 ) 

i i,a=x,y 
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Here i = (i x ,i y ) is the site index, and a is the unit vector 
along a = (x,y) direction, tpi represents the four anni¬ 
hilation operators of an electron on site i. The model 
parameters A, B , C, and D can be experimentally con¬ 
trolled and a is the lattice constant. Specifically, two 
important physical parameters are coupling strength be¬ 
tween inverted bands A and the mass M. a and r 
are Pauli matrices in spin and orbital spaces, respec¬ 
tively. We consider the long-range disorder potential V 

at fi with Vi = EnLiUoeM-K - g| 2 /(2£ 2 )], where 
Uq is uniformly distributed in (—W/2,W/2) with disor¬ 
der strength W and Nj impurities are randomly located 
among N lattices at {ri, r* 2 ,..., rjv 7 },22r— We fix the im¬ 
purity density n = Ni/N = 5% and the disorder range 
£ = 2a, where different n and £ will not significant in¬ 
fluence our results. Since the two spin block are decou¬ 
pled, we only consider the spin-up block in the rest of 


the paper. Because one spin species of QSHI is a QAH 
insulator, our results in this paper are applicable to the 
QAH or Chern Insulator. 

In our numerical calculations, we study the localiza¬ 
tion length Al as well as the dimensionless intrinsic con¬ 
ductance g of the cylindrical sample with width (circum¬ 
ference) L.— which can eliminate the effect of the he¬ 
lical edge states, g is defined as 1/g = 1/gL — 4/N, 
with gL the two-terminal conductance and N the number 
of propagating channelsThe localization length Al 
of the sample is calculated by transfer-matrix method^ 
with the sample’s length 10 6 to 10 7 . The two terminal 
conductance at Fermi energy Ep is evaluated by the 
Landauer-Biittiker formul a 25 ’ 31 ! 32 with a disordered mid¬ 
dle region with size L x L being considered as coupled to 
two clean semi-infinite leads. In addition, we also inves¬ 
tigate the energy level statistics as well as participation 
ratios. For simplicity, the parameters B = 1, C = 0, 
D = 0, and a = 1 are fixed in the rest of paper, where 
we have assumed the particle-hole symmetry (D = 0). 


III. METAL-INSULATOR TRANSITION 


First, we study A = A l/L versus disorder strength 
W by increasing the mass from M = —0.1 to 0.5 at 
a fixed Fermi energy Ep = 0, as shown in Fig. [2] 
These parameters resemble the InAs/GaSb QW k ■ p 
parameters^ Notably, we find the TI-metal-NI transi¬ 
tion in the InAs/GaSb-type BHZ model (unitary system) 
by increasing the mass M. In the NI phase with M < 0, 
A = A l/L decreases monotonously with increasing L in 
FigMa) which indicates all the states are localized. In 
contrast, when M = 0.05 in TI phase (M > 0), the sys¬ 
tem shows one critical (touching) point where A is inde¬ 
pendent of L [see FigJ2](b)], which is consistent with the 
previous studies of a 2D unitary system. ] 16 ’ 18 ’ 19 Due to 
the inverted gap, the TI is robust to weak disorder, and 
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(a) A=0.28, M=-0.1 (b) A=0.28, M=0.05 (c) A=0.28, M=0.38 (d) A=0.28, M=0.5 



FIG. 2: (Color online), (a)-(d) renormalized localization length A = Xl/L versus disorder strength W for different masses M 
and widths L with fixed Fermi energy Ef = 0. The inset of (c) and (d) are the zoom-in of the critical point on left side of 
main panel of (c) and (d), respectively. 


A=0.28, M=0.38, E f =0 




FIG. 3: (Color online), (a) and (c) shows A = Xl/L vs 
disorder strength W near two critical points. The parameters 
are the same as those of Fig. [U)c). (b) and (d) are single 
parameter scaling of A in (a) and (c), respectively, v is the 
critical exponent at critical disorder strength W c . 


the NI phase appears after a certain disorder strength 
W c = 2.3 ±0.13. Therefore, this critical (touching) point 
indicates a direct transition from NI to TI. Surprisingly, 
in Fig. [21(c), when the mass M is increased to 0.38, the 
metallic phase (A increasing with L) appears between 
W c i ~ 1.40 [see the inset of Fig. Efc)] and W c i ~ 7.10. 
This metallic phase is contradict to the common behav¬ 
ior of the 2D unitary class. The unitary system, e.g., 
quantum Hall system, is scaled to localized states except 
at certain critical (touching) points. Moreover, upon fur¬ 
ther increasing M, the metallic phase region becomes 


larger, i.e., for M = 0.5, the metallic phase remains in 
Fig. E](d) between W c \ « 1.40 [see inset of Fig. Hd)] and 
W c2 ~ 8.2. This peculiar metallic phase between the TI 
and NI phases in the InAs/GaSb-type BHZ model is the 
main finding of our paper. 

Next, we analyze the one parameter scaling behav¬ 
ior of the renormalized localization length A = Xl/L 
near the critical points of Fig. He). According to fi¬ 
nite size scaling law, all A are fitted to A(W, L) = 
A c + l a n(W — Wc) n L n t v + boL v near the critical 

point ^ 11 ’ 16 ’ 35 ’ 36 where W c is the disorder strength at 
a critical point, v is critical exponent and y is an ex¬ 
ponent associated with the leading irrelevant operator. 
Here, a n and foo are the fitting parameters. For conve¬ 
nience, we define A '(W,L) = A(W,L) — boL v : with the 
same parameters as A. The best fit is given by minimiz¬ 
ing x 2 statistic y 2 = ~ A(W„, L n )) 2 /cr 2 , where 

N is the number of the data and a n is the error of the 
nth data A n . The finding of the metallic phase and the 
large critical exponent v (almost twice as large as that 
obtained previously) in Figs(3[c) and (d) suggests the 
existence of a possible new universality classi 19 ’ 37- — We 
caution that the standard deviations (x 2 M0 are 3.3 and 
3.5 for Figs[3][b) and (d), respectively, which are some¬ 
what larger than the common value unity and even larger 
size computations are desirable to obtain a high precision 
critical exponent v. 

Up to now, we have found the metallic phase for small 
A = 0.28, which resembles the parameters of InAs/GaSb- 
type BHZ model. To compare with previous studies , 17 ’ 18 
it is necessary to investigate Anderson transition for large 
A cases, i.e., the HgTe/CdTe-type BHZ model. In Figs.Q] 
(a) and (b), where M = 0.38 and Ep = 0, A increases 
from 1 to 2.1, in the region of the HgTe/CdTe-type model 
parameters. The other parameters are the same as those 
of Fig. He). Compared with the case A = 0.28 [see 
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FIG. 4: (Color online), (a)-(d) A = Xl/L versus disorder strength W. The parameters are shown on top of each panel. The 
inset of (c) shows a zoom-in of the critical point on the left side in the main panel. The inset of (d) shows average intrinsic 
conductance (<?) of square sample versus W with same model parameters as (d). 


Fig. He)], the metallic phase disappears for A = 1 and 
A = 2.1 because A hardly changes with the size near 
the touching point [see Figs. H a ) and (b)]. This direct 
transition from TI to NI (i.e. A = 1, 2.1) is consis¬ 
tent with previous study in the HgTe/CdTe-type BHZ 
mode h 17 ’ 18 Meanwhile, the localization length decreases 
rapidly with increasing A from 0.28 to 2.1. For example, 
for a typical unitary case A = 1, the localization length 
A l is only 1^2 times of width of the system L. From 
the above results, it is natural to conclude that the ex¬ 
istence of the metallic phase is highly dependent on the 
magnitude of A. To be specific, the transition Tl-metal- 
NI in the InAs/GaSb-type BHZ model is absent in the 
HgTe/CdTe-type BHZ model. 

Now we consider the influence of Fermi energy Ep 
on the metallic phase. In Fig. He), when A = 0.28, 
M = 0.38, and Ep = 0.08 in the neighborhood of the 
gap center, the transition TI-metal-NI remains almost 
the same as that of Fig. He) with Ep = 0. When the 
Fermi energy is moved outside the band gap E g ~ 0.16, 
i.e., Ep = 0.17, the transition from metal to NI is still 
observed by scaling A and ( g ) in Fig. HH). It is clear 
that ( g) and A increase with the size of the system in 
the metallic phase, they decrease in insulator phase, and 
there is a crossing point (W « 8) at the phase transition. 
To conclude, the localization length scaling and the con¬ 
ductance scaling suggest that the metallic phase could 
survive both inside and outside the band gap. 


IV. ENERGY LEVEL STATISTICS AND 
PARTICIPATION RATIO 

Furthermore, we have verified the existence of the 
metallic phase by studying the energy level statistics and 
evaluating the participation ratio of eigenstates. Ac¬ 
cording to the random matrix theory, the delocalized 
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FIG. 5: (Color online), (a) solid blue line is the variances of 
the energy spacing distributions at different energies Ef with 
the disorder strength W = 4. The red dash line indicates 
the variances of Wigner surmise of unitary ensemble, (b)-(i) 
The histograms of nearest energy level spacings centered at 
different energies Ef = 0-8 with W = 4. The solid red lines 
and green lines represent Wigner surmise of unitary ensemble 
and Poisson distribution, respectively. The model parameters 
are A = 0.28 and M = 0.38 with the sizes 48 x 48. 


and localized states can be characterized by the en¬ 
ergy level statistics! 14 ’ 40 In FigH the histograms of the 
level spacings are drawn at different energies Ep. The 
sizes is 48 x 48 with a periodical boundary condition 
in two directions, i.e., a torus geometry4T 42 In the en¬ 
ergy region Ep = 0-0.2, the histograms [see FigJSjb)- 
(d)] and corresponding variances [see Figj5]( a )] are very 
close to those of the Wigner surmise of a unitary en¬ 
semble p(s) = ( 32 / 7 t 2 )s 2 exp[— (4/7 t)s 2 ]3 Therefore, the 
level correlation is long-ranged which indicates that ex¬ 
tend states may exist. However, near the top of the band 
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( E = 8) [see FigOJi)], the histogram is close to the Pois¬ 
son distribution p{s) = exp(—s), which implies that the 
states are uncorrelated and localized. 

However, in the histograms of the level spacings, it is 
difficult to distinguish the critical point ( or region) from 
the true metallic phase, because they are both close to 
Wigner surmise of the unitary class. To identify the true 
metallic phase, we test the larger sized (LxL = 128 x 128) 
results at fixed energy Ef = 0 with different disor¬ 
der strength W = 1-12. Then we fit the distributions 
to Wigner surmise p(s) = as a exp (bs^) for the whole 
s region and Poisson distribution p(s) = fcexp(Ks) for 
s > 2 , with the fitting parameter a (/3, n) and the 
renormalization parameter a ( 6 , fc)4^ The fitting error 
of a set of data {p(si),p(s 2 ), ...,p(sjv)} is defined as 
<J 2 _ J2? =1 [p(s n )/p(s n )-l} 2 /N, where p(s n ) is the fitting 
value at s„ and N is the number of data. It is noted that 
if energy level spacings data follows the Wigner surmise 
with a = 2 and /3 = 2 , a metallic phase exits; however, 
if the data only follows this distribution in the small s 
region but violates in the large-s region, the critical be¬ 
havior emerges j 43 ' 44 

Figur^Hl shows a logarithmic plot of p(s) for the dis¬ 
tributions of the unfolded level spacings s— In Figsj 6 ][c) 
and (d), for W = 4 and 6 the fitting parameters a and 
/? are very close to the Wigner surmise of the unitary 
ensemble, a = /3 = 2. Therefore, the level correlation is 
long ranged. On the other hand, for W = 1.4 and 8.2 
[see FigHHb) and (e)], although the fitting parameters a 
and /3 are very close to the Wigner surmise of the uni¬ 
tary ensemble, the large s region is not well described 
by the Wigner surmise. Instead, the tail of the large s 
region is clearly fitted to Poisson distribution. The hy¬ 
brid of the Wigner surmise statistics and Poisson statis¬ 
tics behavior at a critical point is coincident with the 
level statistics results at the 3D Anderson metal-insulator 
transition ) 44 In this case, the level correlation is finite- 
ranged, which manifests as a crossing over from (long- 
ranged correlated) Wigner surmise statistics to (uncor¬ 
related) Poisson statistics. For W = 1 and W = 12 [see 
Fig © a) and (f)], the level spacings clearly deviated from 
the Wigner surmise of the unitary ensemble for the whole 
s region and the energy levels will become totally uncor¬ 
related in thermodynamic limit, identifying an insulating 
phase. In conclusion, the level statistics results strongly 
support the existence of a metallic phase between W = 4 
and W = 6 at Ep = 0. This is coincident with the re¬ 
gion indicated by finite size scaling results of localization 
length (see Fig.l). 

The participation ratio characterized the spatial ex¬ 
tension of the eigenstates. Next, we will investi¬ 
gate the participation ratio under different disorder 
strength. The participation ratio is defined as R = 
i \ai\ 4 /{N Y^iLi \ a i\ 2 ) 2 where <Xj is the wave function 
at lattices i and N is the number of the lattice . 31 ’ 46- — R 
has its maximum 1 for one single Bloch wave, reaches a 
finite value typical around 1/3 for a disordered extended 
state, and approaches ~ 1/IV for a localized state. In 


FiglUb) and (c), the (highly possible) extended states 
characterized by the peaks R ~ 0.6 are found near the 
bulk gap at W = 1, and then come into the gap upon 
increasing the disorder strength W = 3. As a result, 
the metallic phase in the gap indeed comes from the ex¬ 
tended states near the gap, which originates from the 
nearly tv Berry phase. At last, the peak of R diminishes 
upon increasing the disorder strength [see FigJ7][d)-(f)] 
which indicates that the states are localized at W = 14 
(R < 1/3). 


V. PHASE DIAGRAM AND DISCUSSION 


We summarize our results into three phase diagrams: 
W-M , W-A, and W-E F , as shown in Fig. [I] In Fig.QJa), 
the metallic phase emerges from TI phase at M ~ 0.1 and 
expands with increasing M from 0.1 to 1. On the other 
hand, in Fig. [D[b), the metallic phase exists for small 
A cases and disappears for A > 0.56. In other words, 
the transition TI-metal-NI can exist in the InAs/GaSb- 
type BHZ model but not in the HgTe/CdTe-type BHZ 
model [see the dash lines in FigOJb)]. Furthermore, in 
Fig.QJc), the metallic phase spreads over both inside and 
outside the gap (E g = 0.16). The phase diagram Fig.[Qc) 
indicates that the metallic phase inside and outside the 
gap have the same physical origin. 

In a 2D system, the metal phase is related to WAL, 
which adds the positive correction to the /3 function and 
leads to such phase in the thermodynamic limit s 24 ’ 49 In 
general, the WAL exists in SOC systems of a symplectic 
ensemble^ and Dirac systems with a tv Berry phase, e.g., 
graphene and helical surface states of 3D TL 50 ’ 51 In the 
present case, the detailed analysis of symmetry classes 
and Berry phases can help us to understand the peculiar 
metallic phase in the inverted band BHZ model. We first 
analyze the symmetry classes of the BHZ model. The 
spin-'I' part of the BHZ Hamiltonian in momentum space 
is H 0 (k) = A{k x r x + k v Ty) + (M — Bk 2 )r z , where A, B , M 
and r have the same meaning as in Eq. ©. When the 
mass term M—Bk 2 F = 0, the system satisfies pseudo-time 
reversal symmetry, i.e., Ho(kp) = t v Hq (—fc_p)r y , which 
is similar to the massless Dirac particles in graphene. 
Therefore, the system is approximate symplectic and 
shows WAL in the region near M — Bk F = 0 j 2 O i 21 On the 
contrary, when |M — Bk 2 F \ » A\k x + ik y |, the two or¬ 
bital bands are nearly decoupled, thus the system belongs 
to orthogonal ensemble approximately and shows weak 
localization ! 20 ! 21 Therefore, the existence of the metallic 
phase is dependent on the system parameters. Moreover, 
the Berry phase of the BHZ model at Fermi energy Ep 
reads : 22 1 26 


7 (Ef) 


(1 


(M - Bk 2 F ) 

y]A 2 k 2 F + {M — Bk 2 F ) 2 


( 3 ) 


where ftp is the momentum at Ep. The Berry phase 7 
monotonously increases with kp rising from zero. 7 = 0 
at fcp = 0, 7 = tv at kp = y/M/B, and 7 — > 2n while 
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FIG. 6: (Color online). Distributions of the unfolded level spacings s under different disorder strength (a) W = 1 to (f) W = 12 
are fitting by Wigner surmise p(s ) = As a exp(Bs^) for whole s region and Poisson distribution p(s) = J 4exp(Bs) for s > 2. 
The model parameters are A = 0.28 and M = 0.38 at Ef = 0. 


A=0.28,M=0.38 
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FIG. 7: (Color online). The participation ratio R as functions 
of energy E at the disorder strength (a) IF = 0 to (f) W = 14 
at the size 48 x 48. 


&F tends to oo. Since 7 can vary from 0 to 2-7T, the sys¬ 
tem can show both WAL and weak localization which 
depends on the parameters. While M/B » A 2 /B 2 , 

k g = \]TT — — \JM/B , where k g is the momen¬ 

tum at the conduction-band bottom and the valence- 


band top, and then the Berry phase 7 ~ n. In this 
case, the system exhibits WAL and the metallic phase 
exists between TI and NI phases (see Fig. [[]). For the 
case of the InAs/GaSb-type BHZ model, i.e., A = 0.28, 
B = 1, and M = 0.38, M/B » A 2 /B 2 is well satisfied, 
the metallic phase consequently appears. In contrast, for 
the HgTe/CdTe type model, i.e., A = 1, B = 1, and 
M = 0.38 with M/B < A 2 /B 2 , k g = 0 which is far away 
from yjM/B. At k g = 0, the Berry phase 7 = 0, which 
leads the weak localization behavior and the disappear¬ 
ance of the metallic phase. 


VI. DISCUSSION AND CONCLUSION 

The Berry phase argument is well applicable to the 
present numerical simulations and it is also in accor¬ 
dance with previous investigations . 20 ’ 21 When the sys¬ 
tem is “approximate symplectic”, electrons interfere as a 
symplectic system on the length scalings smaller than 
Im = vpTM- Here vp is the Fermi velocity, tm = 
^(A 2 kp + (M — Bkp ) 2 )/(M — Bk 2 F ) 2 is the TRS-breaking 
scattering time and r is the elastic scattering time j 20 ’ 21 
In this case, Im acts as large-size cufoff for WAL cor¬ 
rection resembling the dephasing length, and the WAL 
correction dominates if system width L « Im - 20 ' 21 In 
our numerical simulation, for example, when Ep = 0.17 
in Fig.4(d), Im = vfTm = IsTm/f ~ 300, where the 
mean free path l s ~ -^/l/n ~ 4 and tm/t ~ 75. Here we 
have used the impurity density n = 5% and the Fermi 
vector kp ~ 0.605. Since Im » L the system width in 
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the numerical simulations, the TRS-breaking scattering 
time tm is not important and the system resembles sym- 
plectic with metallic behaviors. This metallic behavior 
can show up in mesoscopic systems^ and may apply to 
recent transport experiments^. For larger systems with 
L > Lm, the scaling behavior should be interesting, how¬ 
ever this is beyond our present numerical capability, and 
thus is left for further study. 

In summary, we investigated the Anderson metal- 
insulator transition in QSHI and found different localiza¬ 
tion behaviors depending on model parameters. Notably, 
the transition TI-metal-NI likely exists in InAs/GaSb- 
type systems but not in HgTe/CdTe-type systems. The 
peculiar metallic phase, which originates from the Berry 
phase 7 r near the band gap, contradicts the common view 


of the Anderson transition behavior of the 2D unitary 
class. 
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